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Figure 1. a) bright field LEEM image (field of view: 15 ^m, e^-energy: 2.68 eV) of 
the two graphene phases and uncovered Ir(lll) at room temperature, b) Zoomed out 
PEEM image of the same area as in a) (field of view: 25 /im) The area used for the 
microdiffraction pattern in d) is illuminated by the electron beam (bright circle), c) 
microdiffraction pattern of the parallel graphene phase (electron energy: 45.9 eV) d) 
Zoomed out PEEM image of the same area as in a) (field of view: 25 /im) The bright 
spot marks the area used for the microdiffraction pattern in e) e) microdiffraction 
pattern of the rotated graphene phase (electron energy: 45.9 eV) 
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1. Parallel and rotated graphene phases in LEEM and PEEM 

We observe two domains of graphene on Ir(lll). In [1] a) a bright field LEEM image 
is shown. Graphene in the upper area is imaged brighter than graphene in the lower 
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area, and there is a gulf of bare iridium cutting in from the right. The images b) and 
e) are photo emission electron microscopy (PEEM) images of the same region. Iridium 
is imaged black, graphene from the upper area of a) is imaged darker than graphene 
from the lower area. The bright circle displays the position of the electron beam for 
the LEED images in c) and e). The LEED patterns show that in the case of bright 
graphene in LEEM (dark in PEEM), there is a rotation of graphene's [1120] direction 
with respect to the substrate's [101] direction whereas in the dominating phase - darker 
in LEEM (bright in PEEM) the dense packed directions of graphene and Ir(lll) are 
parallel. 



2. Delamination around heptagon-pentagon pairs 




Figure 2. STM topograph of graphene on Ir(lll) (20 nm x 20 nm). The moire rows 
around the bulge are marked with bright lines. 

The STM topograph in [2] shows a bulge in the graphene layer. The extra ending 

line of moire maxima (three on one side of the protrusion, but four on the other) indicate 

an extra row of atoms terminating in a pair of heptagon-pentagon carbon rings [1]. We 

speculate that the bulge is a nucleus of a wrinkle originating from such a point defect. 
[1] J. Coraux, A. T. N'Diaye, C. Busse, and T. Michely. Structural coherency of graphene on Ir(lll). 
Nano Lett., 8:565-570, 2008. 
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Abstract. Upon cooling, branched line defects develop in epitaxial graphene grown 
at high temperature on Pt(lll) and Ir(lll). Using atomically resolved scanning 
tunneling microscopy we demonstrate that these defects are wrinkles in the graphene 
layer, i.e. stripes of partially delaminated graphene. With low energy electron 
microscopy (LEEM) we investigate the wrinkling phenomenon in situ. Upon 
temperature cycling we observe hysteresis in the appearance and disappearance of 
the wrinkles. Simultaneously with wrinkle formation a change in bright field imaging 
intensity of adjacent areas and a shift in the moire spot positions for micro diffraction 
of such areas takes place. The stress relieved by wrinkle formation results from the 
mismatch in thermal expansion coefficients of graphene and the substrate. A simple 
one-dimensional model taking into account the energies related to strain, delamination 
and bending of graphene is in qualitative agreement with our observations. 



1. Introduction 

The new material graphene receives currently an enormous attention for its exciting 
properties. At the heart of the scientific interest are the consequences of graphene's 
unique band structure arising from its lattice symmetry and its monoatomic thickness 
[1]. The high mobility of electrons in graphene and the strong electric field effect 
foster work to realize graphene based electronics [2]. Moreover, the use of graphene 
for conducting transparent electrodes [3l H], to realize photosensitive transistors [5], 
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ultracapacitors [Hj , or a new class of catalytic and magnetic materials through templated 
cluster growth has been suggested [71 [8] . 

Although the exciting electronic properties of graphene have been explored 
mainly by transport measurements for devices built on flakes of exfoliated graphene 
on Si02, there appears to be consensus that for future scientific exploration and 
technological applications epitaxial growth of high quality graphene over large areas 
[9l[IOl[ni|12lll3l[Ta[l5]isa prerequisite. 

For any technological application it is of utmost importance to avoid or at least 
control the defects in graphene associated with the epitaxial growth process. Such 
defects may result from growth obstacles caused by substrate steps [ini IH] or by the 
coalescence of finite sized graphene domains [13] . 

Branched line defects are present in mono- or multilayers of continuous graphene at 
room temperature after high temperature (>1000K) epitaxial growth on several metals 
and on SiC. Their nature has been under debate. While some authors attributed the 
branched line defects to carbon nanotube formation [161 [13 [H [IE] , current research 
evidences that they are wrinkles in the graphene layers on SiC [201 1211 [221 [1] as 
well as on metals [231 [2^ [25] . For wrinkles in graphene on SiC very recent atomic 
resolution scanning tunneling microscopy (STM) data provide unambiguous evidence 
for the wrinkle interpretation [22] . Very recent LEEM and STM data provide evidence 
for the presence of wrinkles also in monolayer graphene on Ir(lll) [25]. Wrinkle features 
have also been observed in layered transition metal dichalcogenides [26]. Among other 
suggestions several authors attribute the occurrence of wrinkles to the difference in the 
thermal expansion coefficients of graphene and its support [231 [201 [2ll [221 [25] . Consistent 
with that mismatch is the observation of compressive strain in epitaxial graphene at 
room temperature [271 [28l [7] . This compression was found to vary on a length scale of 
less then 300 nm [29]. Although in situ investigations could provide deeper insight into 
the wrinkling phenomenon, such investigations are missing till now. 

Here we not only support the wrinkle interpretation of the branched line defects 
by atomically resolved scanning tunneling microscopy (STM) imaging of wrinkles on 
metals, but also gain a detailed insight into wrinkle formation by in situ LEEM imaging 
and micro diffraction. Most noteworthy, LEEM and micro diffraction "see" not only 
wrinkle formation but also the structural and electronic changes in adjacent /im-sized 
areas within the graphene. A model is developed, which qualitatively agrees with our 
observations. 

2. Methods 

We examined epitaxial graphene on Ir(lll) and Pt(lll). Graphene has been grown 
epitaxially by chemical vapor deposition of ethene (C2H4) at elevated temperatures in 
ultra high vacuum. Scanning tunneling microscopy (STM) was carried out at room 
temperature, low energy electron microscopy (LEEM) imaging was done at variable 
temperature. Growth and imaging was performed in ultra high vacuum without any 
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transfer outside the vacuum. 

3. Results and discussion 

3.1. Wrinkle formation 

Figure [T]a) shows a bright field LEEM image of a graphene fiake on Ir(lll). The fiake 
has a diameter of ~ 6 /im. Branched line defects on the fiake which develop upon cooling 
to room temperature form a network of dark lines, much darker than the substrate step 
structure which can faintly be seen in the image as well [31]. The STM topograph 
in figure [T]b) shows graphene islands on Pt(lll). The branched line defects are also 
present. A typical line defect is shown in figure [T]c). It crosses the image diagonally 
and diverges in two at the bottom of the image. It roughly follows the direction of the 
two monoatomic steps of the underlying Ir(lll) substrate. The defect is about 3nm 
high, and thus much higher than a substrate step as shown in figure [T]e). Its width is 
a few nanometers as well. On low line defects as shown in figure [Ud), profile in figure 
[Hf), it is possible to achieve atomic resolution on the ridge. The atomic rows over the 
defect are continuous. The dense packed rows cross the wrinkle roughly perpendicular. 
The fact that these structures are present only on the graphene fiakes and never on the 
uncovered part of the surface corroborates the assumption that they are indeed wrinkles 
and not nanotubes on the sample, as has been proposed previously [Ml [IH [TBI [19] . The 
continuity of the atomic rows also indicates that the elongated structures are not formed 
at ruptures where the islands edges roll or bend up. 

The LEEM images in figure [2]a) show epitaxial graphene on Ir(lll) fully covering 
the field of view. As visible in micro diffraction (see supplement) the darker graphene 
domain on the left is rotated with respect to the substrate [25]. The left image [figure 
[2]a) I] has been taken at high temperature (HOOK), close to the growth temperature 
(T^row=1320 K) of graphene. No wrinkles are observed. During cooling wrinkles appear 
and spread all over the field of view as visible in figure [2]a) II. Upon reannealing close to 
the growth temperature, the wrinkles disappear again [figure [2]a) III]. Faint dark lines 
due to steps are present at all temperatures. The time and temperatures the images 
have been recorded at are marked with green dots in the temperature vs. time diagram 
in figure [2]b). 

The appearance of a wrinkle is an instant process within the time resolution of our 
measurement of 1 s, while the decay of the wrinkles is a gradual process. Wrinkles decay 
at slightly higher temperatures than they form. The red squares in figure [2]b) indicate 
the formation and decay temperatures. This also refiects in a hysteresis of the average 
lattice parameter of graphene as measured by spot profile analysis LEED measurements 
[32]. 

Graphene has been grown at several temperatures and the onset temperature of 
wrinkle formation has been recorded. The onset of wrinkle formation is measured as 
the occurrence of the first wrinkle in a field of view of 10 fim. Soon after the occurrence of 
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Figure 1. a) LEEM image (field of view: 10 /im, electron energy: 2.8 eV) of a 
graphene flake on Ir(lll) at room temperature, b) STM topograph of graphene flakes 
on Pt(lll) (3 /i7Ti X 3 /i7Ti). The image is differentiated and appears as if illmninated 
from the left. The inset shows the same image, with the graphene flakes highlighted in 
red. c) STM topograph (240 nm x 240 nm) of a full layer of graphene on Ir(lll) with 
a wrinkle. The bright line corresponds to the profile given in e). d) STM topograph 
(7nm X 7nm) of a wrinkle of low height in graphene on Pt(lll) in atomic resolution. 
The bright line corresponds to the profile given in f). e) Profile of the topography of 
the wrinkle in graphene on Ir(lll) shown in c). f) Profile of the topography of the 
wrinkle in graphene on Pt(lll) shown in d). STM images have been processed using 
the WSxM software [3D]. 
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Figure 2. a) Three bright field LEEM images (field of view: 10 /im, electron energy: 
2.7 eV) of Ir(lll) fully covered with two differently oriented domains of epitaxial 
graphene. Images I and III have been taken at 1110 K while II has been imaged 
at 410 K. b) Temperatm'e evolution measured at the sample with a pyrometer (black 
curve), and at the sample holder with a thermocouple (grey curve). The green dots 
mark the points where a) I-III were recorded. The red squares mark the onset of wrinkle 
formation upon cooling and the disappearance of the last wrinkle upon heating, c) 
The difference AT of the onset of wrinkle formation and the graphene growth 
temperature Tgrow based on the pyrometer measurement is plotted as a function of 



the first wrinkle additional wrinkles emerge and continue to form even until the sample is 
cooled down to room temperature. The first wrinkles appear after a cooldown of 410 K 
± 40K regardless of the growth temperature [see figure [2]c)]. Assuming an average 
difference of the thermal expansion coefficients of Ir and graphene of 7.15 x 10~^ a 
temperature difference translates to a mismatch in expansion. During a cooldown by 
410 K iridium shrinks by 0.33% while a graphene layer only shrinks by 0.03% [331 IM] . 
The remaining 0.3% have to be taken up by compression or wrinkling of the graphene. 
This suggests a stress driven process for wrinkle formation. 
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Figure 3. a) Micro diffraction pattern of grapliene on Ir(lll) near the specular beam 
with first order moire spots, b) Relative length ki/ki dm'ing cooldown (see text), c) 
Effect of strain relief in reciprocal space. Only the component of the reciprocal lattice 
vectors ki which is parallel (light colors) to direction of strain relief (thick black line) 
is influenced, d) Intensity of the specular beam during cooldown. 



3.2. Lattice expansion 

Direct evidence for abrupt strain relief during cooldown is shown in figure [31 Panel a) 
shows a micro diffraction pattern near the specular beam at 19.9 eV and resulting from 
an area of about 1.4 fim diameter. At this energy, three of the six moire spots are bright. 
The reciprocal lattice vectors of the moire fcj are the difference of the corresponding 
reciprocal lattice vectors of graphene and Iridium. For small relative changes, their 
behavior reflects the relative length changes and angular changes in reciprocal lattice of 
the graphene layer amplified by factor of ten [35] . The evolution of the relative distances 
ki/ki of the bright moire spots to the specular beam is displayed by red, green, and blue 
lines in figure [3] b) |36] . Here, ki denotes the time averaged value of ki during the first 
25 s, i.e. prior to the sudden decrease at about t = 30 s. This decrease amounts to 
si = (3.6 ± 0.4)%, S2 = (2.0 ± 0.4)%, and S3 = (1.2 ± 0.5)%, where Si = 1 - k/k'^ 
with k^ being the average for the length of the reciprocal moire lattice vector between 
t = 32 s and t = 42 s. We attribute this decrease to a transition from the fiat state into 
the wrinkled state (see section 3.3). The wrinkled graphene is relaxed and has a larger 
lattice parameter and thus shorter reciprocal lattice vectors ki [37] . The three ki shrink 
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by different amounts. This is well understandable given the wrinkle is a linear defect 
which can only relax stress in one direction. Assuming all three ki to be of equal length 
and making angles of exactly 120° prior to wrinkling, we can estimate the direction 
of stress relaxation. Each vector fcj can be split into a component which is parallel 
to the direction of strain relief and one which is perpendicular to this direction. Only 
the parallel component is affected by the strain relief. Let (pi be the angle between 
the lattice vector ki and the direction of strain relief, ki and k'- the averaged lengths of 
the reciprocal lattice vectors before and after the transition, fcj y = fcjcos^j the average 
component parallel and fc, ^ = ki sin 0j the average component of ki perpendicular to 
the direction of strain relaxation and prior to the relaxation and finally c = k'- \\/ki^\\ the 
reciprocal space compression upon wrinkling. This leads to 

1 - = (1 - c^) cosV. 

for all three reciprocal linearly dependent lattice vectors. With the measured values for 
^ this equation can be numerically solved. This leads to a factor c = 96% corresponding 
to a strain relief of 0.4% in the graphene layer. The angle between the direction of 
compression and ki is roughly (pi = 10° as illustrated in figure [3] c). This implies that 
the wrinkling took place roughly perpendicular to the dense packed rows of the graphene 
layer consistent with what is expected from STM data as the dominant orientation of 
wrinkles (cf. figure [T]d) or [22]) The intensity of the specular beam abruptly increases 
simultaneously with the relaxation of the lattice as shown in figure [3] d). 

3.3. Local stress evolution 

The abrupt intensity change of the specular beem upon wrinkling allows us to monitor 
the local extension of stress relaxation through LEEM imaging. Figure Hla) shows LEEM 
images of graphene on Ir(lll). The sample is partially covered by graphene prepared at 
1110 K and the sample has been cooled down to 560 K within one hour. In the course of 
this cooling, some wrinkles have already formed, especially in the right hand side of the 
field of view. The images I and II have been measured subsequently with a delay of 1 s. 
They capture a single event of wrinkle formation on the graphene patch near the left 
border of the image. Simultaneously with wrinkle formation the brightness increases 
in the affected area. This change is visualized in figure |4]a) III which is the difference 
between the images in II and I. It shows that the formation of a wrinkle does not only act 
nanoscopically at the line of delamination but it rather has an impact on a mesoscopic 
scale (in this case 4/im^). 

The change in intensity integrated over the regions marked in figure Hla) I is shown 
in figure |4]b). This change takes place at different times, locations and to different 
extents. We interpret the increase in image intensity as an effect of the relaxation of 
the graphene lattice. This provides an explanation for the observation of locally varying 
compression in epitaxial graphene [29] . 
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Figure 4. a) Two subsequent LEEM images of graphene on Ir(lll) recorded during 
cooling (field of view: 10 ^m, electron energy: 2.5 eV) and the difference of their 
intensities with enhanced contrast. Between (I) and (II) a wrinkle forms on the 
peninsula on the left. This is shown enlarged and contrast enhanced in the insets. 
The difference image shows that the intensity has increased locally in the course of 
wrinkle formation, b) The intensity integrated over the regions marked by colored 
boxes in a) I over time, with linear background subtraction. Whenever a wrinkle is 
created, the brightness increases abruptly. The roman numbering indicates where a) I 
and a) II have been recorded, c) LEEM I(V)-curve. Blue: imwrinklcd graphene, red: 
wrinkled graphene, purple: difference (x20). 

To discuss the origin of the (0, 0)-spot intensity increase upon wrinkhng we consider 
figurelUc). It shows I(V)-curves of the (0, 0) spot of flat and wrinkled graphene on Ir(lll) 
and their difference. The curves were taken from one sequence of images recorded at 
intermediate temperatures where wrinkles have just started to form and some parts 
of the surface are still unwrinkled. The difference curve has two maxima. Both are 
correlated with maxima in total intensity resulting from constructive interference of 
electron waves reflected from graphene and the Ir surface. For a distance of 3.4 A [38] 
we expect maxima at 3.2 eV and 13 eV. Therefore we speculate that the intensity changes 
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are due to structural changes, such as a change in the spacing between the graphene layer 
and the Ir substrate triggered by the relaxation of the graphene lattice. Evidently these 
structural relaxations are accompanied by changes in the graphene electronic structure 
which additionally may affect the I(V)-curve. 

During a wrinkling event many atoms are displaced. On Ir(lll) and Pt(lll) 
graphene forms an incommensurate superstructure [35l [39] . That implies a low barrier 
for sliding of graphene on the surface, because for every atom, which loses energy by 
moving out of its optimum binding configuration another atom gains energy. The small 
flakes grown by temperature programmed growth [13] barely show wrinkles. A graphene 
flake smaller than the average separation of two wrinkles just expands if the compression 
gets large enough to overcome the barrier for sliding. For smaller islands, this barrier 
may be larger, due to edge effects becoming more relevant. Accordingly there is residual 
strain in such flakes |35j . 

Wrinkle patterns from repeated cooling and heating cycles at one sample spot are 
similar. This suggests that wrinkles nucleate at preexisting features. We find in our STM 
data spots of delamination which are centered at heptagon-pentagon pairs of carbon 
rings (see supplement) |10]. This seems reasonable as these defects induce additional 
local stress into the graphene lattice. We thus speculate that heptagon-pentagon defects 
are sites of wrinkle nucleation. 

3.4. Model 

The wrinkle formation can be described in a one-dimensional continuum model. When 
the substrate and graphene cool down, the graphene has to compensate for the thermal 
misfit AL/L resulting from the difference of the thermal expansion coefficients. Either 
compression ALc/L or the formation of a wrinkle AL^/L can compensate for that misfit 
{AL = ALc + ALyj). 

To calculate the energy of compression we use the separation between two 
wrinkles L = 260 nm as estimated from experiment. With an atom density of 
n = 36.2 atoms/nm^ and an elastic modulus of F = 56 eV per atom [H] the compression 
energy per nm of wrinkle length can be expressed as Ec = ^{AL^/ L^L ■ Y ■ n. The 
energy necessary for wrinkle formation consists of two contributions: first, there is the 
bending of the graphene layer. An estimation for this contribution is available from the 
study of single walled carbon nanotubes [H], giving the bending energy per atom in a 
nanotube of radius r as = a/r^ + b/r^ with the empirical parameters a = 0.99eV/A^ 
and b = 12.3 eV/A^. Second, bonds between graphene and the substrate are stretched 
or broken, where the graphene flake delaminates. Since the estimated height of graphene 
on Ir(lll) is comparable to the interlayer distance of graphite, we use the exfoliation 
energy of graphite as an assessment for the binding energy of an atom in a graphene 
layer on Ir(lll) in this case. It has the value of Cbfl = 0.052eV [l2l H3] The strength 
of the van-der-Waals binding energy between a particle and a surface decreases with 
the cube of the distance, so we calculate the change in binding energy of an atom in an 
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Figure 5. a) A map of the energy per nm of a wrinkle ii^tot with respect to 
uncompressed flat graphene according to a one-dimensional continuum model. The 
shape of the wrinkle has been modeled as four arcs of circles (see inset). The lower 
horizontal axis shows the misfit of graphene and Ir the upper axis shows the 
corresponding temperature difference AT. The vertical axis represents the fraction of 
the misfit which is compensated by wrinkling ^jf- instead of compression, b) The 
graph shows the energy of a wrinkle per nm for a mismatch of 0.7% (thick black line) 
and the contributions it consists of. There is a local minimum for the flat configuration 
where all the energy is stored in the form of compression E^, while in the global energy 
minimum most of the energy is stored in a wrinkle as bending and reduced bonding, 
c) Energy barrier E^.w foi' wrinkle formation. The system can gain energy by wrinkle 
formation (i?gain) with respect to compressed flat graphene, if the misfit is larger than 
0.65%. 
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elevated part of the graphene sheet with the height of z instead of zq above the substrate 
as Cb = Cbfi (1 - z^/z^). 

With and Eb as the integral of 6^ and Cb over the atoms in a nm of wrinkle the 
resulting energy cost for a nm of wrinkle compared to fiat relaxed graphene sums up to 



and leads to a complex variation problem for the shape of the wrinkle. 

Here we assume a simple shape for the wrinkle, consisting of four equal arcs of a 
circle with radius r and an opening angle a as shown in the inset of figure [5]a). The 
model contains r (or a) as a free parameter which is optimized for minimum energy Efot 
for each combination of AL/L and AL^/L. 

In figure [5]a) i?tot is plotted for this model as a function of thermal misfit AL/L 
(lower horizontal axis) and fraction of strain accommodated in a wrinkle AL^j/AL 
(vertical axis). Moving along the horizontal axis from left to right corresponds to cooling 
of the sample. The according temperature difference is indicated on the upper horizontal 
axis. In the lower part of the diagram most of the misfit is taken up by compression 
while in the upper part the misfit is compensated for by a wrinkle. For misfits below 
0.3% there is exactly one optimum configuration: The graphene layer is compressed 
and there is no wrinkle. As the misfit increases (temperature decreases), a second local 
minimum in energy emerges. Nevertheless, the unwrinkled compressed state still is 
favorable. For AL/L > 0.65% of misfit, a situation rendering about 80% of the misfit 
subject to wrinkle formation, is optimal. 

A cut through the map at constant misfit AL/L = 0.7% is shown in figure Ob). 
There is a local minimum for the fiat configuration where all the energy is stored in the 
form of compression Ec, but the optimum configuration is the formation of a wrinkle, 
which contains most of the energy in the form of reduced bonding to the substrate (Eb) 
and bending of the graphene {E^)- Figure [5]c) illustrates the relationship of the two 
energy minima and the barrier in between. The gray line shows the difference in energy 
of the unwrinkled state and the wrinkled state -Egam- For a misfit below 0.47% there is no 
minimum for the wrinkled state, above 0.47%, there is a local minimum, but its energy 
is higher than that of the uncompressed flat state. Only for compressions above 0.65%, 
when the gray line enters the negative region, the system can gain energy by forming 
a wrinkle. Still, there is an energy barrier to overcome, which allows the system to be 
trapped in the local minimum explaining the sudden and abrupt formation of wrinkles. 
This is consistent with the hysteresis for wrinkle appearance and disappearance. 

Although our one-dimensional model explains all qualitative features observed, the 
model prediction overestimates the experimentally observed critical misfit for wrinkling 
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formation nearly by a factor of two. Certainly a full two dimensional analysis may lead 
to somewhat different numbers - wrinkle formation is likely to be eased by biaxially 
compressed graphene. Also the wrinkle separation L, the binding energy Ehfi and our 
simple model for the wrinkle shape carry significant uncertainties. 

As wrinkles are large scale defects, it would be desirable to suppress their formation. 
One way to achieve this could be to reduce the amount of total thermal misfit by growing 
graphene at the lowest possible temperature and inserting an intermediate annealing 
step to remove the defects prior to cooldown. Combining temperature programmed 
growth and chemical vapor deposition it appears possible to achieve high quality 
graphene at a growth temperature of only 1000 K [Hj. Also grazing incidence keV 
ion erosion removing exclusively protruding wrinkles followed by annealing could lead 
to continuous graphene with less or no wrinkles. A third approach could be to increase 
the energy for bending. This could be accomplished by evaporating at high substrate 
temperatures a thin film with a thermal expansion coefficient similar to the substrate 
on top of the graphene. This cover layer would have to be bent as well, for the graphene 
to form a wrinkle. 

4. Conclusion 

In conclusion, we demonstrated how LEEM can be used to monitor strain relaxation in 
situ. It was possible to develop a consistent picture of wrinkle formation on graphene 
linking wrinkle formation with inhomogeneous residual strain. The development 
of wrinkles appears to be a serious problem for all methods of growth of weakly 
bound epitaxial graphene, as all require high temperatures. We hope the improved 
understanding of wrinkle formation achieved here will contribute to a future solution of 
this problem. 
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